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The impact of coherence on the nonlinear optical response and stationary transport is stud- 
ied in quantum cascade laser structures. Nonequilibrium effects such as pump-probe signals, the 
spatio-temporally resolved electron density evolution, and the subband population dynamics (Rabi 
flopping) as well as the stationary current characteristics are investigated within a microscopic 
density-matrix approach. Focusing on the stationary current and the recently observed gain oscil- 
lations, it is found that the inclusion of coherence leads to observable coherent effects in opposite 
parameter regimes regarding the relation between the level broadening and the tunnel coupling 
across the main injection barrier. This shows that coherence plays a complementary role in sta- 
tionary transport and nonUnear optical dynamics in the sense that it leads to measurable effects in 
opposite regimes. For this reason, a fully coherent consideration of such nonequilibrium structures 
is necessary to describe the combined optical and transport properties. 

PACS numbers: 78.20.Bh,74.78.Fk,42.65.-k,73.63.-b 



I. INTRODUCTION 

Quantum cascade lasers (QCLs)i are semiconduc- 
tor heterostructure lasers employing the transitions be- 
tween quantized intersubband levels in quantum well 
structurcs^i^i^i^i^ and act as a source of radiation in the 
terahertz (THz) /mid- infrared regime. The laser consists 
of multi-quantum well periods comprising the electron in- 
jector and the optically active region. These periods arc 
repeated tens or even hundreds of times over the length 
of the structureii To drive the electrons through the sam- 
ple, an external bias is applied. Scattering processes and 
optical recombination between the conduction subbands 
in the doped structure within as well as between periods 
lead to stationary electronic occupations out of equilib- 
rium. While on the technological side, this light source 
can be used for spectroscopy in the fields of environmen- 
tal detection or medicinejii^ it offers on a fundamental 
ground an interesting model system to study intersub- 
band charge dynamics in a structure where the optical 
and the transport properties arc closely interrelated. 

The first semiconductor heterostructure laser operat- 
ing in the intersubband regime was realized by Faist et 
al. in 1994ii Since then, many types of QCLs of dif- 
ferent design have been built and optimized, see, e.g., 
Refs.00 for an overview. The QCL has been the subject 
of extensive theoretical research. The stationary proper- 
ties of QCLs were studied by a rate equatio n^":^^ and 
a Boltzmann-type approac h^^i^'^i^^i^^i^^i^^ as well as a 
quantum theory employing both nonequilibrium Green's 
functionsi^ii^i2£i2ii^ and density-matrix theoryi^^i2iiS 
Here, the gain, the current-voltage characteristics, and 
the stationary charge distributions have been established. 
First results regarding the nonlinear optical properties 
such as optically induced subband population dynam- 
ics have been presented by us within a density-matrix 



theoryi^ 

One central result of these studies is that coherence 
can play an important role in the determination of the 
stationary current, requiring a fully microscopic the- 
ory of the current including non-diagonal density-matrix 
elements i^°i^^'^^ Here, the often applied rate and Boltz- 
mann equation approaches, based on Wannier-Stark hop- 
ping (WSH), fail. Experimentally, indications of coherent 
charge transport have been obvscrvcd in the oscillatory 
gain recovery in pump-probe experiments of mid-infrared 
QCLs2S>2£i^i2S as well as recently in THz structures.— 
This oscillatory behavior has been attributed to reso- 
nant tunneling through the injection barrier of the laser. 
These observations suggest that coherent effects might 
also become visible in optics in the time regime be- 
yond the light-matter interaction. However, different 
studicsMi^I reveal a simple relaxation in the gain recov- 
ery, showing that not all samples exhibit this coherent 
effect. 

Motivated by these investigations, this paper is focused 
on the role and importance of coherence in the interplay 
between ultrafast optical dynamics and stationary trans- 
port in quantum cascade structures. We investigate the 
regimes where coherence is of relevance in the combined 
optics-transport system. To this end, we present a fully 
microscopic theory describing the dephasing and tunnel- 
ing processes in a QCL structure. In order to systemat- 
ically investigate the optical and transport regimes, we 
focus on two quantities describing the coherent and the 
incoherent evolution within the system: the tunnel cou- 
pling 20, between the two states across the main injection 
barrier and the level broadening F of the states. We sys- 
tematically vary the width of the main tunneling barrier 
in the QCL structure, thus establishing a relation be- 
tween these two central quantities. Using this relation, 
we investigate the importance of coherence in the calcu- 
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lated signals: In the transport regime, we compare a rate 
equation (WSH) approach with a fuUy microscopic cur- 
rent theory to investigate the regime of vaUdity of both 
models. In the optical regime, we consider pump-probe 
calculations, complemented by the spatio-temporally re- 
solved electron density evolution and the subband popu- 
lation dynamics due to strong ultrafast excitation. Using 
these nonlinear results, aspects of the nature of the op- 
tically induced charge dynamics and the importance of 
coherence can be addressed 

Comparing the optical and transport results, we find 
that the influence of coherence is observable in oppo- 
site parameter regimes: the inclusion of coherence in 
stationary transport calculations becomes important for 
2n < r, where the coherence between subbands limits 
the current flowing through the structure, while in the 
nonlinear optical gain dynamics, the inclusion of coher- 
ence drives the observed oscillations in the gain recovery 
for 20 > r. Combining the two results, we thus find 
that it is necessary to consider a fully coherent theory in 
order to understand the combined optics-transport prob- 
lem independent of the parameter range of operation. 



II. THEORY 

The QCL heterostructure is modeled as a multi- 
conduction subband system, each period comprising an 
injector and an active region. Within the effective mass 
and envelope function approximations^ the wave func- 
tions, assumed to be separable and infinitely extended in 
the quantum well plane, are given by 

where A is the in-plane quantization area, x and k are 
the in-plane spatial and momentum vector, respectively, 
S.i{z) is the envelope function in growth direction, and 
Uj:j^g(x) the Bloch function taken at the band edge. 
There are several natural possibilities to choose the wave 
functions S,i{z) to describe the system. Since we consider 
the optical response of the system as well as the WSH 
model in the current calculation, we choose the Wannier- 
Stark (WS) basis which diagonalizes the heterostructure 
potential as well as the externally applied bias. This 
offers an intuitive physical interpretation of the optical 
transitions and the scattering between single-particle WS 
states as well as of the electron counting between approx- 
imate system eigenstates. It should be noted that, even 
though the physical observables are in principle indepen- 
dent of the choice of basis, different results are expected 
for different basis choices due to the necessary approxi- 
mations applied in many-particle problems. In the fol- 
lowing, the theory is applied to the THz QCL discussed 
in Ref . [33. The band structure of this QCL in the original 
design is shown in Fig. [1] with the important subbands 
marked. This structure is considered throughout this pa- 
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FIG. 1: (Color online) Band structure of the THz QCL from 
Ref. for the injection barrier width b = 5.5 nm under res- 
onance condition. The subbands involved in the laser transi- 
tion (4,1') as well as the injector subbands (5,2') are marked. 



per, considering different injection barrier widths b. The 
parameters used in the calculations are found in Table [D 
In order to describe the ultrafast nonlinear optical and 
stationary transport properties of the laser, it is prefer- 
able to consider a fully microscopic treatment of the scat- 
tering mechanisms in the structure. In this work, we 
consider the interaction of the electronic system with lon- 
gitudinal optical (LO) phonons via the polar (Frohlich) 
coupling as well as with the ionized doping centers in the 
laser. Even though Coulomb scattering can be treated 
in the same way^^^ we do not consider this interaction 
here. For the relatively lowly doped THz laser consid- 
ered here, it was shown that scattering with impurities 
typically dominates over Coulomb scattering4S Acous- 
tic phonons are excluded since they induce only very 
long scattering times in quantum wells and act mainly as 
a low-energy intrasubband thermalization mechanism4i 
Interface roughness scattering is excluded due to the difR- 
culty to quantify it microscopically, while alloy scattering 
is small in structures with a binary material for the wells 
(GaAs) and therefore neglected4^ The included scatter- 
ing mechanisms are treated in a density-matrix corre- 
lation expansion approach within a second order Born- 
Markov approximation42iii We do not include renormal- 
izations of the subband energies due to scattering in the 
form of principal values (see Sec. Ill B 31 for a discussion 
of this approximation) ; it is assumed that these are small 
and approximately constant for the different transitions, 
leading only to an absolute energy shift without physical 
consequences (see Ref. for a discussion of scattering- 
induced energy renormalizations in quantum wells). 



A. Hamiltonian 

To derive the dynamical equations, we divide the 
Hamiltonian of the system into three parts: 

H — Hq + Hcl-light + Hscatt- (2) 



3 



TABLE I: Structural and material parameters used in the 
calculations. 



material system 

doping density n2d 

well electron mass m* 

LO phonon energy fio^LO 

high-frequency permittivity £00 

static permittivity Sa 



GaAs/Gao.isAlo.ssAs 
2.945x10^° cm-2 
0.067 mo 
36.7 meV 
10.9 
12.9 



The first part, 

Ho = HsL + H, + Ho.ph = '^tk'^lk'^.k + ^q^q^q' 

tk q 

(3) 

describes the kinetics of the electrons in the hetcrostruc- 
turc potential Hgi^ as well as the externally appUed bias 
i/j and the kinetics of the phonons -ffo.ph- The electronic 
part is diagonal due to the choice of basis, a^- (a -) 

ik ik 

denotes the creation (annihilation) operator of an elec- 
tron in subband i vi^ith quasi-momentuni k and energy 
^ik ~ £i + fi-^fc^/(2"T.i) and b]^ (b^) the creation (annihila- 
tion) operator of an LO phonon with three-dimensional 
quasi- momentum q and energy fkuq = huLO- Subband 
nonparabolicity is neglected, and the subband effective 
masses are assumed to be constant, nii = m* , where m* 
is the effective mass of the well material. 

Hei-Mghi. describes the interaction of the system with 
a coherent classical light source. The polarization of the 
light field is chosen in the direction of the dipole moment, 
i.e. in the growth (z-) direction, and the field is assumed 
to be a spatially homogeneous Gaussian pulse, 

A(t) = A{t)e,,A{t) = Ao exp{-(t/T)2} cos(tJLt), (4) 

with the laser frequency ljl and Gaussian pulse dura- 
tion T. The pulse area is defined via the envelope of the 
pulse4^ Under the assumption of a homogeneous excita- 
tion, the Hamiltonian is momentum-diagonal and reads 



i?ei-iight = eA{t) Y J2 



ik jk 



(5) 



with the elementary charge e > 0, the coupling ele- 
ments = ^{i\[p^ /m{z)+eA{t) / {2m{z))]+Ii.c.\j) , and 
the momentum operator pz = {h/i)dz with the space- 
dependent effective mass m{z). 

Finally, iJscatt describes the scattering processes: 
electron-LO phonon interaction i/ci-ph as well as scat- 
tering with ionized doping centers i/ci-imp- 



1. Electron-LO phonon interaction 



by 



The electron-LO phonon coupling Hamiltonian is given 



'J fe,q 



where q is the in-plane projection of q. The coupling 
matrix element is given by the Frohlich coupling 



q 

9u = 



2eoV \e 



1 



1 



.1/2 



e(q) • q 



(7) 

Here, Eg and Soo are the static and high-frequency permit- 
tivity, respectively, V is the quantization volume, e(q) is 
the displacement unit vector, and q±^ the projection of 
the momentum q in growth direction. 



2. Interaction with ionized doping centers 

Since the quantum cascade structure is doped, it is 
necessary to take into account the interaction of the elec- 
trons with the ionized doping centers. Typically, either a 
barrier or a well in the QCL is doped, making it necessary 
to distribute the ions in this layer. Here, we treat the in- 
teraction following Ref. [l^, where the dopant density is 
distributed on several 5-sheets located at zi in the doped 
barrier /well. The ionized doping centers are considered 
as classical scattering centers. The Hamiltonian for the 
interaction is given by 



i^eW.p = EEE^^('^ 
'J k,q I 



at a . (8) 

ik+q jk ^ ' 



with the screened electron-impurity interaction potential 



Here, A is the screening constant and I labels the position 
of the randomly distributed individual ion (5-distributions 
in the growth direction. For the screening, we use the 
static limit of the Lindhard formula for a homogeneous 
electron gas 41 



B. Dynamical equations 

1. Equation structure and approximations 

We derive the dynamical equations using a correla- 
tion expansion within a density-matrix approach in sec- 
ond order Born-Markov approximation, applying a bath 
approximation for the LO phonons424i This leads to 
equations of motion for the microscopic polarizations 
f - 

ij.k 



ik 



a^-a.-)(i 7^ j) and the subband occupations 
On the level of the phonon- (impurity )- 



ik jk 



ik 



assisted coherences (a'^^ba^a ((a^^Vl,„((f)a - )), re- 

* ik ^ jk" ik jk" ' 

spectively, we neglect the interaction with the opti- 
cal field since these terms are of higher order in the 
coupling!^ We assume the system to be homogeneous in 
the plane perpendicular to the growth direction z; then. 
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we can restrict to a density matrix diagonal in the in- 
plane momentum: (aj^a^.^,) = f.^^khk'- 

The general dynamical equations are found in the ap- 
pendix. While for the interaction with the doping cen- 
ters, the complete equation structure is considered, only 
the terms linear in the density matrix are taken along 
for the electron-phonon interaction due to numerical rea- 
sons. This approximation is justified if we assume that 
we are working in the regime of nondegenerate electron 
gases, where n^j ^ 1 which is typically fulfilled in these 
QCL structures. 

We should stress at this point that, due to the spatial 
extension of the Wannier-Stark wave functions, it is es- 
sential to consider the whole set of matrix elements in 
the calculations. Testbed calculations considered only 
certain sets of matrix elements, e.g., we restricted to 
the diagonal/nondiagonal scattering terms gf^gf*, con- 
sidered typically in nonbiased (equilibrium) quantum 
well systems,— i^i^i^ or alternatively to terms where the 
overlap of the wave functions is the main argument, i.e. 
terms such as gf^gf* along with the above set of terms. 
Both versions lead to nonstable results in some param- 
eter regimes where the coherences are important. It is 
thus essential to consider the full set of matrix elements 
when considering the properties of the QCL at resonance 
condition and for large injection barrier widths, i.e. small 
tunnel coupling across the injection barrier. 

When considering the WSH model, the relaxation pro- 
cess are restricted to Boltzmann scattering given by 

^(0) ^ _pout,^ _^ pin _ ^ Y /^QN 
ik tk ik ' ik^ ik' ^ ' 

The in- and out-scattering rates F'" and r°-* are 

ik ik 

functions of the subband occupations n^p and are given 
for both the electron-phonon and electron-impurity 
interaction in the appendix. 

The system temperature enters the calculations 
through the scattering rates which include the phonon 
distribution, assumed to be given by an equilibrium Bose- 
Einstein distribution riq = [exp(?kjq/fcsT) — 1]^^, and 
the screening parameter A for the interaction of the elec- 
trons with the doping centers, see also Ref. HO- Here, the 
approximation is made that the lattice and the electronic 
temperature are the same. 



2. Periodic boundary conditions 

In order to describe the extension of the periodically 
coupled structure, it is necessary to apply appropriate 
boundary conditions in the growth direction. Since the 
inclusion of coherence between states of different periods 
is important to describe the nonlinear optical dynam- 
ics due to the spatial extension of the wave functions,?*' 
we consider a nearest neighbor approach where the co- 
herence and wave function overlap between two adjacent 



periods (n, n') is taken into account. We checked this ap- 
proximation by showing that |/„^„, ^1 , gZ.j„, , V;'„^„, (q) ~ 
for \n^n'\ > 2, where {gl^„,,Vi^„, (q)) denotes 

the coherence (coupling element) between the states \i) 
and \j) in the periods n and n', respectively. We apply 
translational invariance of the density matrix between 
periods)^ 

finjn'M ^ /i(ri+l)j(n' + l),fe' "in,fc ~ "'i(ri-(-l),fc- (-'"-'-) 

For the wave functions and the band edge energies Si , we 
apply a coordinate shift and a bias drop, respectively, 

£,in{z) = £,i{n+l){z + ipcr)i ^in = ^ 6£'-Lpor, 

(12) 

where Lpcr is the period length and E the applied electric 
field. In the following, the index i = (i,n) is taken as a 
composite index. 

3. Shortcomings of the model 

Even taking along the whole set of equations within the 
approximations discussed above, we partially encounter 
negative occupations (see Fig.[2|) which is known from 
the treatment of Redfield-type equations such as consid- 
ered hercj^i^i^ However, this typically does not lead to 
unphysical results if we focus on averaged observables. 
For all dynamical calculations, we obtain physical values 
for the current density, the total electron densities in each 
subband, and the gain spectrum. However, the broaden- 
ing in Eq. (|13p may become negative for regimes where 
coherences play an important role due to the negative oc- 
cupations. This may be due to the fact that (i) the scat- 
tering rates as such are not observables and (ii) the simple 
formula for the broadening used here, which is typically 
applied in the literature, does not take into account sta- 
tionary coherence. In order to guarantee strictly positive 
values for the averaged scattering rates to calculate the 
broadening in Eq. (|13p . which we can use as a measure 
of the lifetime of the Wannier-Stark states, we perform a 
Gaussian smoothening of the stationary distributions and 
use these in the calculations of the broadening F. Fig- 
ure [2] shows an example of this smoothening for two dif- 
ferent barrier widths: While for small barrier widths and 
low temperatures, the negative occupations are negligi- 
ble, and thus the corresponding averaged scattering rates 
are always positive, the smoothening becomes important 
for large barrier widths where the negative occupations 
can lead to negative values. Thus, the values obtained for 
the broadening F should be viewed as approximative val- 
ues. Since we are concerned with qualitative results only, 
this is justified. For the cases where the negative occu- 
pations are negligible, and thus the averaged scattering 
rates strictly positive, the smoothening induces changes 
of the broadening of 2 — 3%. We would like to stress 
again, however, that all dynamical calculations and all 
results except for Fig. [5] are performed without artificial 
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k^/(2 m) (meV) ^ k^/(2 m) (meV) 

FIG. 2: (Color online) Gaussian smoothening of the subband 
occupations n^k for the two main subbands (injector 5, upper 
laser state 1) at T = 10 K: (left) h = 4.5 nm, (right) 6 = 11.0 
nm. 



modifications of the calculated data, such as any kind of 
smoothening. In addition, it should be noted that the 
calculations here are performed at the point of "maximal 
coherence" with respect to the tunnel coupling across the 
injection barrier, i.e. at resonance (cf. Sec. IIII A]) . Away 
from resonance, the influence of coherence, and the cor- 
respondingly related problems such as non-positivity, are 
strongly reduced. 

The reason for the quite large negative values for the 
occupations partially encountered may be explained by 
the fact that the system under consideration is a very 
complicated one: for the nearest neighbor coupling con- 
sidered in the calculations, we have 35 independent coher- 
ences, and thus there is a very strong interplay between 
the different frequencies of the coherences which can de- 
stroy the strict positivity of the occupations. For small 
barrier widths and low temperatures, the general dynam- 
ical equations (|A.ip and \K.h\ reduce to the Boltzmann 
dynamics given in Eq. (jlOp . Here, any negative occu- 
pations observed in our full calculations are negligible 
(«ik > -10~®) [see, e.g.. Fig. [^left)]. In the case of the 
Boltzmann scattering dynamics, i.e. including only the 
diagonal elements of the density matrix, the stationary 
occupations remain strictly positive for all barrier widths 
and temperatures. 

The specific oscillation energies of the peaks can be 
explained by the energetic structure of the system: the 
peaks correspond to different combinations of the bias 
drop per period eFLpcr and the LO phonon energy 
E'lo- In Fig- Elright), e.g., the large peaks correspond 
to the energy eFLpcr — i?LO = 56.40 meV —36.7 meV ~ 
19.7 meV, while the smaller peak positions are given by 
2eFLpor — Si^LO = 2.7 meV. The latter are weaker, since 
they represent a higher order process of iterative trans- 
port and scattering. For the two subbands, the peaks are 
shifted by AE w 0.1 meV, corresponding to the tunnel 
splitting between the two states, Ey ~ E^ w 0.1 meV. 

As stated above, there is a strong interplay between 
different coherences of sharp energy. For T = 50 K and 
a certain barrier width regime, we need to incorporate 
additional damping in order to obtain convergence. In 
the calculations, we do this by adding a phcnomcno- 
logical, low-energy intrasubband scattering mechanism 
with an energy which is uncommensurable with the LO 



phonon energy {E « 0.5 meV) and use the LO phonon 
coupling element with an increased strength. We should 
note that this damping is of a purely phenomenologi- 
cal character and not a physical scattering mechanism 
as implemented here. The concerned barrier widths are 
marked in Figs. O and (5] The problem with convergence 
may be a consequence of the combination of neglecting 
low-energy scattering channels and applying the Born- 
Markov approximation within the density-matrix theory. 
Another explanation of this convergence problem is that 
scattering-induced energy renormalizations are neglected 
at this order of the perturbation expansion. This is 
an approximation often applied within the Born-Markov 
approximationj ^^i"^^ Since the renormalizations, like the 
scattering rates, are temperature-dependent, this could 
explain why the convergence problems only appear for 
certain temperatures. On the other hand, since the reso- 
nances depend on the specific energy structure, and thus 
on the barrier width &, problems are expected to appear 
only for certain barrier widths, as witnessed in the re- 
sults. 

Including an artificial damping as discussed above, we 
find that the value of the stationary current is approx- 
imately independent of the scattering strength. This 
has already been observed in earlier studies of the QCL 
(Refs.Hllil) and attributed to the fact that a low-energy 
dissipative scattering mechanisms such as LA phonons 
leads mainly to a thermalization within the subbands, 
but has no strong effect on the actual value of the station- 
ary current. We adopt the validity of this statement here 
and take the obtained values of the current as physical re- 
sults. On the other hand, it is clear that the inclusion of 
additional large damping leads to a strong overestimation 
of the total broadening as calculated in Eq. (flSl) . Since 
the relation between the broadening F and the tunnel 
coupling 2f2 is used to discriminate qualitatively differ- 
ent regimes, the inclusion of additional scattering mech- 
anisms for the concerned barrier widths (which typically 
leads to a larger broadening) docs not lead to qualita- 
tively new regimes of interest (see Fig. [S]). It is for this 
reason that one can safely disregard these values in the 
discussion of the broadening and the tunnel coupling as 
is done in the following. 



III. NUMERICAL RESULTS 

In this section, we apply the theory of Sec. [ll] to the 
THz laser from Rcf. [s^. This structures has been investi- 
gated with respect to its transport in the stationary laser 
regime as well as its gain properties'^ In the regime of 
ultrafast nonlinear optics, Rabi oscillations were consid- 
ered recently)^ 

First, the systematic setup is discussed, along with the 
calculation of the central quantities used for the further 
discussion of the tranport and optical properties. The 
stationary nonequilibrium due to the different scattering 
mechanisms is determined. This is necessary for the dc- 
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FIG. 3: (Color online) (left) Injector (5) and upper laser state 
(1') showing an anti-crossing at the main tunneling barrier at 
exact resonance of the corresponding Wannier states. The 
lower laser state (4) is also shown, (right) Tunnel splitting 
energy AE — {E^ — E-^i \ for varying bias and different barrier 
widths h determining the resonance bias at the center of the 
anti-crossing and the tunnel coupling 2Q. ~ AE. 



termination of the stationary value of the current, but 
also as a starting point for the consideration of the opti- 
cal response. Then, the WSH approach is compared with 
a fuU current calculation to investigate the importance 
of the inclusion of coherence in the stationary transport. 
Finally, the nonlinear optical response of the structure 
is considered, focusing on pump-probe signals as well as 
the spatio-temporally resolved electron density and the 
optically induced population dynamics, to study the role 
of coherence in the optical regime. 



A. Systematic setup and stationary state 

In order to systematically investigate the importance 
of coherence, we vary the width of the main tunneling 
barrier connecting the injector and the active region of 
the laser. When resonant tunneling between the injec- 
tor and the upper laser state is small, the WS states are 
approximately localized in the injector and in the active 
region. For a certain applied bias, the WS states become 
delocalized across the injection barrier and an anticross- 
ing occurs, where the injector and the upper laser level 
form a pair of binding/ antibinding states. The system 
is assumed to be in resonance at the center of the anti- 
crossing, i.e. for the bias at which the splitting energy 
between the two respective levels has a minimum. Here, 
the level splitting AE equals twice the tunnel coupling 
2f2 between the localized injector and the upper laser 
state, see Fig. [31 

The QCL as a multi-subband nonequilibrium system is 
a very complex structure, and thus a straight-forward an- 
alytical determination of the stationary population and 
coherence distributions due to scattering is not possible. 
We thus numerically determine the initial conditions to 
be stationary solutions for f^^ ^ and n^j: without the op- 
tical field by solving all scattering contributions from an 
arbitrary state of fixed total population until a steady 
state is reached. Figure [3] shows the process of deter- 
mining the initial distributions by starting from subband 
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FIG. 4: (Color online) Population relaxation dynamics us- 
ing the full set of dynamical equations for a barrier width 
of 6 = 5.5 nm at T=10 K, 50 K, determining the stationary 
nonequilibrium of the system. Inset: Short-time evolution of 
the subband populations. 



populations Ui = 2/ A'^^n^^ given by single-subband 
Fermi distributions with equal populations in each sub- 
band. An approximate stationarity is reached on a time 
scale of a few picoseconds, while long time intraband 
redistributions can last up to several nanoseconds (not 
shown). For small barrier widths, i.e. large tunnel cou- 
pling f2, the relaxation is described well by a Boltzmann 
relaxation rt*"?. For larger barrier widths, the influence 

ik 

of the coherence on the stationary populations cannot 
be neglected any longer and becomes important which is 
discussed in detail in the next section. 

Having determined the stationary distributions, we 
can calculate the scattering induced broadening of the 
states. The level broadening of the considered states 
is approximately given by the mean of the two tunnel- 
split Wannier-Stark states, F = (Fi' + F5)/2, where the 
broadening is determined by the stationary Boltzmann 
out-scattering rate. 



A 



ik ik 



,(0) 



(13) 



for the stationary occupations At these low tem- 

peratures, the broadening is typically dominated by the 
impurity scattering, which is strongly temperature de- 
pendent due to the screening length^ and thus increases 
strongly from 10 K to 50 K, while the interaction with 
phonons is essentially given by spontaneous emission. 

Figure [5] shows the calculated values of and F as 
a function of the barrier width h. As is expected, the 
broadening remains roughly constant for varying barrier 
widths, since the scattering rates are not strongly influ- 
enced by the width, while the tunnel coupling across the 
barrier approaches zero for b oo. For T = 10 K, the 
two energies are roughly the same for 6 ~ 6.0. This is 
the regime which is physically interesting, and thus bar- 
rier widths slightly larger and less than this value are the 
focus in the following investigations. For T = 50 K, the 
broadening of the states is much larger then the tunnel 
coupling for all considered barrier widths, and thus no 
transition between different regimes is expected. 
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FIG. 5: (Color online) Tunnel coupling 2Q, and mean broad- 
ening r = (Fs + Fi/)/2 of the injector and the upper laser 
state for different barrier widths 6 and at T = 10 K, 50 K." 

B. Stationary results: Current calculations 

The stationary distributions far from equilibrium de- 
termined in Sec. IIII Al are the starting point for the op- 
tical dynamics discussed later in this chapter. Here, we 
consider the accompanying evolution to the steady state 
value of the current flowing through the structure. In or- 
der to investigate the influence of coherence on the sys- 
tem, we focus on two different approaches to the cur- 
rent calculation: the first restricts to the occupations 
n-jt to consider a rate equation approach - the so-called 
Wannier-Stark hopping model - and the second takes the 
coherences f. . r into account. 

We start with the rate equation approach, where in 
effect the electrons are counted as they cross an interface 
at a fixed point in the growth direction of the structure 
due to scattering between states which arc localized at 
different spatial parts of the structure. The WSH current 
is just given by the application of Fermi's golden rule to 
each set of states. For the electron-phonon interaction, 
it is determined via (ji-q = riq) 

^^^^^^ = EE - ± ^LO) 

X ("q + ^ T - ^w>.fe(l - ^fc+,-). (14) 
where denotes the absorption/emission of a phonon, 



^ For the barrier widths marked by dashed lines (open circle in 
the current density), the inclusion of an additional low-energy 
scattering mechanism is necessary to achieve convergence of the 
current (see the discussion in the text). The value of the current 
density obtained in this way is approximately independent of the 
strength of the additional coupling when determined from the 
transient, in accordance with earlier results™!^ The resulting, 
unphysically large total broadening F does not lead to qualita- 
tively new regimes of interest and thus not considered further. 



da = —ezii is the expectation value of z for the WS state 
and the sum z,j is carried out over all states. The 
factor 2 arises due to spin degeneracy and L ~ NLpa- is 
the length of the structure with TV periods of length Lpor. 
The current is fully determined by the diagonal elements 
of the density matrix n.^, considering jumping of elec- 
trons from one state to the next, and reaches its steady 
state value as determined by the relaxation in Fig. [1] 
when restricting to Boltzmann dynamics. An analogous 
expression is found for the elastic impurity scattering. 
While Eq. pi]) can be motivated by hopping of electrons 
between different positions za and Zjj, it is in fact an 
approximation for the full current driven by coherences, 
Eq. HID, see Ref.H. 

When deriving the current microscopically from 
the current operator J = —e/{2■mo)^^x,t)[p + 
eA{t)]^{x, <) + H.C. , one arrives at another picture of the 
current, where it is the coherence between states which 
carries the current flowing through the structure. Here, 
the current is given completely via the nondiagonal ele- 
ments of the density matrix f^j ^, which in growth direc- 
tion yields 

JfunW = -2^^^«.,4.£, (15) 

»i k 

where Vij = ^{i\[pz/'m{z) + eA{t)/m{z)] + H.c.|j) are 
the velocity matrix elements and the factor 2 again 
arises due to spin degeneracy. For the calculation, 
the polarizations f^j g are initially set to zero. The 
finite occupations in the subbands lead to a stationary 
current which is fully determined by the scattering in 
the structure. 

We now carry out current calculations using both the 
WSH and the full approach for varying barrier widths. 
The result is shown in Fig. [51 The WSH current remains 
approximately constant for varying barrier widths. This 
is due to the fact that at resonance, the injector and 
the upper laser level form a pair of binding/ antibinding 
states whose nature hardly changes with b. Thus, the 
hopping current is not affected by the barrier width, see 
also Ref. [13. The slight increase of the current may be 
attributed to the impact of the other levels which change 
slighly with decreasing bias, which is needed to keep the 
two tunnel- coupled states in resonance. 

In contrast, the full calculation [Eq. based on the 
coherences shows the expected drop of the current with 
barrier thickness. We first focus on the case T = 10 K. 
For small barrier widths (where 2il > F), the two sta- 
tionary current results are approximately the same and 
show a qualitatively similar behavior, i.e. the slope of the 
two curves are almost parallel. In this regime, the WSH 
and the full model both describe the stationary current 
well. The slightly smaller value for the full model stems 
from the fact that only the main injection barrier is in- 
vestigated systematically here. Resonant tunneling at 
the other barriers in the structure is not considered, so 
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is shifted to lower barrier widths. It should be noted 
that for all barrier widths here, the relation 20 ^ F is 
fulfilled. 

It is thus found that in the stationary current calcu- 
lations, the inclusion of coherence becomes important if 
20 < r under resonance. For 20 > F, the WSH ap- 
proach to the current calculation is a good approxima- 
tion, and thus the inclusion of coherence is negligible. 
For 20 < F, only the full model shows the physically 
expected decrease of the current as the tunnel coupling 
decreases since the coherence between the injector and 
the upper laser state becomes important to describe the 
localization of the eigenstates. The decrease of the tun- 
nel coupling across the main tunneling barrier limits the 
current flowing through the structure. This limitation 
is not reproduced by the WSH approach. It should be 
noted that this result has already been discussed before 
in Ref. [13 for the QCL and in Ref. [H for superlattices. 



C. Ultrafast optical dynamics 



FIG. 6: (Color online) Stationary current for both the WSH 
and the full calculation at resonance for varying barrier widths 
at T = 10 K, 50 K." 



that a similar argumentation as applied in the following 
can also be applied to them. Thus, it is expected that 
a constant offset can occur, probably due to the extrac- 
tion barrier where the two delocalized states (lower laser 
level, extraction level) are also close to resonance. 

For large barrier widths (where 20 < F), the full cur- 
rent decreases until it almost vanishes in the limit of large 
barrier widths. This decrease is the physically expected 
behavior, as the growing injection barrier width restricts 
the current flowing across the barrier and thus through 
the whole structure. For larger barrier widths, the grow- 
ing localization of the charge in the injection region leads 
to the fact that the Wannier-Stark states which are de- 
localized across the injection barrier no longer constitute 
a "good" basis to describe the states of the system. For 
this reason, the coherence between these two states be- 
comes important in the stationary state, whereas it is 
negligible for smaller barrier widths. This can be seen in 
the relaxation dynamics as well, where for large barrier 
widths, the influence of the coherences on the popula- 
tions, mainly of the injector and the upper laser state, 
becomes very important, and thus a pure Boltzmann re- 
laxation description fails for the stationary populations. 

For the case of T = 50 K. the situation is different. For 
small barrier widths, the drop of the full current is more 
pronounced than for T — 10 K, showing that even for 
those barrier widths the WSH model for the current fails. 
We should note that this difference is not as pronounced 
as would be expected; however, since we are interested 
in qualitative results, it suffices to find that the onset 
of deviation between the WSH and the full calculation 



We now turn to the nonlinear optical response due to 
an ultrafast external perturbation of the laser structure 
from stationarity. Here, the situation is not as clear as in 
the transport case. Obviously, the light-matter interac- 
tion as such is a coherent process as the vector potential 
couples to the nondiagonal elements of the density ma- 
trix f^^ j:. However, it is more interesting to consider the 
dynamics after the passage of the pulse, where the signal 
is typically dominated by incoherent scattering. 

As discussed in Sec. IIII Bl coherences between 
Wannier-Stark states are important in the regime of large 
barrier widths and high temperatures, where the station- 
ary states are coherent superpositions of the WS states; 
they describe the tunneling through the barriers and thus 
the return to the stationary state. This is expected to 
remain important in the optical dynamics in order to 
describe the return to the steady state after optical exci- 
tation. However, the comparatively large scattering with 
respect to the tunneling (see Fig. [5]) is expected to de- 
stroy all observable coherent effects in this regime. The 
question arises whether coherent optical effects after ul- 
trafast optical excitation can be observed in the opposite 
regime, where coherences were shown not to be impor- 
tant in the stationary transport. Recently, optical pump- 
probe measurements in mid-infrared^'2£'^i2^ and THz 
QCLs2^ have shown an oscillatory behavior attributed 
to (coherent) resonant tunneling through the injection 
barrier. In this section, this effect is investigated with 
respect to the relation between the tunnel coupling and 
the level broadening. 

We begin by considering the linear response of the 
system to characterize the laser. Then, we analyze the 
nonlinear optical response in form of the experimentally 
studied pump-probe signals to study the regime of impor- 
tance and the role coherence plays in the optical signal. 
We complement this study by looking at the optically 
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FIG. 7; Dependence of the absorption spectrum on the injec- 
tion barrier width at T = 10 K. 



induced charge as well as the subband population dy- 
namics. 

The linear absorption of the structure is calculated via 
the complex susceptibility 



1 5J{uj) 



(16) 



eq ijj'^A{uj) ' 

SJ{t) = J{t) — Jq denotes the change of the current den- 
sity induced by the optical field with vector potential 
A{oj), where Jo is the stationary value of the current den- 
sity determined in the previous section and J{t) is given 
by Eq. (jl5p . From this, the absorption a{uj) ~ wlmx((x>) 
is calculated. 

In Fig. [71 the absorption spectrum of the gain transi- 
tion is shown for different injection barrier widths. Due 
to the anti-crossing at resonance, the spectrum shows a 
double-peak structure which can be resolved for small 
barrier widths where the tunnel splitting is sufficiently 
large. For larger barrier widths, the two peaks merge to 
form a single resonance for sufficiently small tunnel split- 
ting energies compared to the dephasing of the transition. 
Due to this splitting, it is a priori unclear which of the 
transitions is to be considered the optical gain transition. 
In the following considerations, the system is excited res- 
onantly on the transition yielding the larger peak gain 
which here is the lower gain transition frequency, and 
the accompanying optical dynamics is investigated. The 
second peak appearing for b = 9.0 nm at AE « 16 meV 
arises from an enhanced dipolc moment between the two 
resonant states l',5 and the extractor state 3. 

When increasing the barrier width further, the coher- 
ence between the injector and the upper laser state be- 
comes important, so that the stationary state of the sys- 
tem is in a linear superposition of the two WS states, 
localized in the injector region. This coherence can lead 
to a strong decrease or even a vanishing of the gain. To 
illustrate this, it is helpful to consider a simple three-level 
system consisting of the two resonant states (injector 5 
and upper laser state 1') and the lower laser state 4. We 
assume that the coherences between the lower laser state 
and other two states is small compared to the station- 
ary tunneling coherence at the time of the probe pulse 



to, which is taken as a (5-pulse, i.e. |/4i'(to)|, 1/45(^0)! ^ 
l/i'5(*o)l- Solving the semiconductor Bloch equations, 
the gain at the laser frequency is then given by (under 
the assumption |i?5 — Ey] |i?5 — £^4]) 



41'(»1' 



ru) + 2di5divRc [/i'5(to)]} 

with the constant level densities rii . If the system is now 
in the state \'^) |5) + |1') (o?45 « —d^i/) which is lo- 
calized in the injection region, as is the case for large 
barrier widths, we have 715 « nii . The gain is then given 
by the stationary inversion of the states 715 — 714 as well as 
a further term — 2Re [/i'5(to)] containing the coherence 
between the two resonant states. 



-dls {2(715 - 7^4) - 2Re [/i'5(M]}- 



(18) 



Typically, the first term containing the inversion de- 
termines the optical spectrum. In the presence of a 
strong stationary coherence /i'5(to)i the gain can be 
strongly reduced or even vanish, even in the presence 
of a strong inversion, which is the case for all barrier 
widths considered. The stationary state as a super- 
position of the two WS states is fully localized in the 
injector, leading to a very small dipole overlap with the 
lower laser state which is witnessed in the vanishing gain. 

For the nonlinear optical dynamics, we focus on the 
pump-probe signals which have been experimentally 
measured recentlyi^>2£i^i22iSi2ii^ An ultrafast nonlin- 
ear pump pulse resonant on the laser transition excites 
the sample, leading to a gain saturation at the laser en- 
ergy, and a subsequent weak probe pulse tests the laser 
transition as a function of the delay time between the 
two pulses (see Fig. [51). Corresponding to the experi- 
ments, a 170 fs Itt pump pulse (P) and a (in compari- 
son to all optical transitions of interest) spectrally broad 
probe pulse (Pr, with T = 50 fs) are used. In contrast to 
a standard differential transmission spectrum where the 
complete energy range is recorded;^ the focus is on a 
monochromatic measurement, yielding the absorption at 
a fixed energy for different delay times. The experimen- 
tal separation of pump and probe pulses which is done 
via directional filtering of the response is carried out here 
in such a way that the response to only the pump pulse 
5Jp is subtracted from the response to both the pump 
and probe pulses (5Pp-|-pr, yielding the purely linear re- 
sponse after the influence of the pump pulse whose direct 
response is filtered out. Thus, we consider for the absorp- 
tion 



a{td,uj) ^ Im 



SJp+Pritd,^) — 5Jp{uj) 



UjApy{td,Uj) 



(19) 



with the delay time td, the fixed pump and probe en- 
ergy (D (which is equal to the laser energy here), and 
the probe field Api-. 5Ji[uj) is the Fourier transform 
of the corresponding simulated time signal 5Ji (t) taken 
at the fixed frequency ui. In Ref. [sgI . it is shown that 
the nonlinear absorption according to Eq. describes 
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FIG. 8: (Color online) Pump-probe signals after ultrafast 
strong pumping and subsequent linear probing at the laser 
energy for different delay times td and barrier widths at T = 
10 K, 50 K. Both pump and probe pulses are chosen to be 
resonant on the laser transition. Insets: Pump-probe signals 
for small delay times. 



the signal contribution measured in the probe direction. 
While the signal contains all coherent effects also in the 
limit of delay times shorter than or comparable with 
the dephasing times, additional coherent effects are in- 
cluded for short delay times which are not measured 
experimentally due to directional filtering; these effects 
are expected to be small. Our result reproduces well 
known studies in the probe direction, for instance Ref . [stI. 
The differential pump-probe response is then given by 

exp{-[a(td)-ao(t<i)]i}-l ~ ~a{td) + aoitd) = PP(td), 
where ao(td) denotes the linear response without the 
pump pulse and it is assumed that aL ^ 1 where L 
is the length of the structure. For reasons of presen- 
tation, we plot the negative differential pump-probe re- 
sponse. Then, positive signals correspond to a decreased 
gain compared to the stationary reference gain value. 

The pump-probe response for different barrier widths 
and at T = 10 K, 50 K is shown in Fig. [5] (see also Fig.[n] 
for the corresponding population dynamics caused by the 
pump pulse for b = 4.5, 7.0 nm at T = 10 K). The ultra- 
fast pump pulse saturates the gain transition, leading to 
a strong absorption at this energy (positive pump-probe 
signal). The following return to the steady state shows 
two main features: (i) A decay of the signal on a time 
scale of a few picoseconds. This feature is due to the in- 
coherent scattering which leads to a return to the steady 
state, yielding the original gain for the laser transition, 
(ii) For smaller barrier widths, an oscillatory modulation 
of the signal whose amplitude decreases for growing bar- 
rier widths, while its period of the order of picoseconds 
increases correspondingly. In addition, fast oscillations 
with a period on the order of 0.1 ps are visible which are 



addressed briefly below, but which arc not the focus of 
the article. 

The oscillatory feature is a coherent effect in the dy- 
namics after the passage of the pulse: The pump pulse 
depopulates the upper laser state locally in the excited 
region. At resonance condition, such a locality is repre- 
sented by a coherent superposition of the levels 1' and 5. 
This causes coherent charge oscillations between the two 
superpositions which are localized in the injector and in 
the active region;^ leading to a modulation of the coher- 
ence between the two laser states and thus to an oscilla- 
tion of the gain signal. The oscillation hardly causes any 
changes in the WS subband populations and no oscillat- 
ing inversion as will be discussed later in this section. 
Thus, the oscillation is a coherent effect. The charge 
oscillations are visible for the case that the tunneling 
period Tqsc ^ 1/AE is less than the scattering period 
Tscatt ~ l/F, i.e. Tosc < Tscatt- Let us first focus on T = 
10 K. Here, gain oscillations are strongly visible for b ~ 
4.5 nm, 5.5 nm where 2fl > F, i.e. Tqsc ^ Tscatt- Already 
for b = 6.5 nm, the oscillations become hardly visible 
and fully vanish for d > 7.0 nm, where Tosc 2> Tgcatt- 
For very large barrier widths, it is not possible to invert 
the gain transition due to the strong coherence between 
the injector and the upper laser state. The two states 
are now in a superposition state localized in the injec- 
tor, and thus a dipole interaction via the electric field 
between the lower laser state and the upper state is not 
possible since the charge is localized outside the active 
region. Thus, no charge inversion occurs (compare the 
discussion of the gain for large barrier widths in Fig. [7]). 
At T = 50 K, 2n < F, i.e. Tosc > T^catt is fulfilled for all 
barrier widths (cf. Fig. [5]), and thus only very weak gain 
oscillations are found for b = 4.5 nm. Again, for large 
barrier widths, no inversion occurs, leading to an almost 
vanishing pump-probe signal. It should be noted that 
the fast oscillations on a time scale of around 300 fs seen 
in the insets of the pump-probe signals for small delay 
times id ^ 1 ps can be attributed to tunneling due to 
the coherence between the lower laser state and the two 
resonant states which show an energy splitting of AE « 
13 meV, corresponding to an oscillation period of Tqsc ~ 
300 fs. Due to the changing applied bias and the corre- 
spondingly changing laser energy, the oscillation period 
changes slightly between Tqsc ~ 200-300 fs. 

Experimentally, gain oscillations have been observed 
in a mid-infrared laser in Refs. and in a THz 

laser in Ref.lH. In the case of the mid- infrared laser, pro- 
nounced gain oscillations were found up to relatively high 
temperatures, including a gain overshoot when probing 
close to the resonance. This effect depends strongly on 
the strength of the scattering mechanisms, specifically on 
the depletion of the lower laser subband and the strong 
coupling through the injection barrier, and has not been 
observed in our calculations. According In Ref. [3l|, the 
lifetime of the superposition of the two resonant states 
Tscatt ~ 1 ps which is much longer than the observed os- 
cillation of Tosc ~ 500 fs, and thus 2il > T which is the 
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FIG. 9: (Color online) (left) Spectral function and (right) the 
imaginary part of the lesser Green's function for two different 
barrier widths at T = 10 K. 

regime where gain oscillations are expected. In Ref. IsBI . 
time-resolved pump-probe differential transmission mea- 
surements did not show signs of gain oscillations. As the 
authors argue in the paper, this is due to the very short 
scattering rates due to Coulomb scattering (Tscatt < 100 
fs) which is much smaller than the tunneling oscillation 
period (Tqsc = 517 fs), and thus no gain oscillations are 
expected. 

The interpretation of the oscillations in the pump- 
probe signals being caused by coherent charge transfer 
across the injection barrier is strengthened by considering 
the spatio-energetic structure of the stationary state ob- 
tained from our Green's function modeli^ In Fig.[HlJleft), 
the spectral function A{E, k = 0, z) which describes the 
density of states for vanishing in-plane momentum is 
shown. The peaks roughly correspond to the WS wave 
functions |^i(z)P with an energetic width F^. For b = 
4.5 nm, the binding and antibinding combination of the 
injector and the upper laser level can be resolved around 
-E w 36 mcV, while this is not the case for b ~ 7.0 nm as 
AE < F. The right figures show the lesser Green's func- 
tion Im [G^ {E, k = 0, z)] describing the corresponding 
carrier density. For b = 4.5 nm, the carrier density qual- 
itatively follows the spectral function corresponding to 
two occupied levels in the stationary state. Thus, when 
perturbing the stationarity on an ultrashort time scale, 
a superposition of the two states leads to an oscillatory 
modulation of the gain. For b = 7.0 nm, the carrier 
density does not follow the density of states; only the su- 
perposition of the two states, localized in the injector, is 
occupied at stationarity, and thus no superposition can 
be excited and no oscillations occur. The stationary dis- 
tribution is already a superposition of the two WS states, 
which is also the reason why the inclusion of coherence is 
important to correctly describe the stationary transport 
of the system (see Fig. [H]). 

In addition to the spatio energetic resolution of the sta- 
tionary charge density, it is insightful to consider its spa- 
tiotemporal evolution during and especially after the cx- 
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FIG. 10: (Color online) Spatio-temporal evolution of the elec- 
tron density of the QCL after ultrafast nonlinear optical exci- 
tation at T = 10 K. The pump pulse parameters are the same 
as for the pump-probe results. 

citation with an ultrafast laser pulse. To do this, we 
consider the spatio-temporally resolved electron density 
which is given by 

■y k 

where the time dependence is included in the density ma- 
trix f^j j:{t). In order to simplify the representation and 
to focus on the optically induced dynamics, the difference 
between n{z, t) and the stationary electron density 710(2;), 
An(z,t) = n{z,t) — 71,0(2:), is considered in the follow- 
ing. Similar calculations have been presented in Ref. [23l . 
where however no detailed modeling of the pump pulse 
was done. 

In Fig. [ini the electron density evolution is shown for 
the same parameters of the pump pulse as used for the 
calculation of the pump-probe signal discussed above. 
Shortly after the pulse has passed at i « 0.3 ps, the 
density is depleted in the active region due to the fast 
nonradiative extraction from the lower laser level to the 
extractor state and into the injector. 

After this, two qualitatively different behaviors are 
seen for the two barrier widths. For b = 4.5 nm, the 
resonant tunneling across the injection barrier leads to a 
fast refilling of the upper laser state and thus to an in- 
creased density in the active region at i w 1 ps. At t ~ 2 
ps, the density in the active region again decreases. Thus, 
the oscillation of charge which occurs on a time scale of 
1.8 ps as detected by the probe pulse, see Fig. [5{left), 
can be directly seen in the electron density evolution. 
When considering the electron density without the po- 
larizations, i.e. restricting to the diagonal elements of 
the density matrix for the calculation, or even just with- 
out including the coherences between different periods 
(not shown), this oscillation is not observed at all. 

For the case of 6 = 7.0 nm, a roughly monotonous re- 
turn of the density depletion to the steady state is found, 
and thus no gain oscillations are expected, as verified 
in the pump-probe signals. Due to the weaker tunnel 
coupling between the injector and the active region, the 
charge transfer is slow so that the oscillation period in- 
duced by the tunnel coupling Tqsc is larger than the scat- 
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FIG. 11: (Color online) Optically induced population dynam- 
ics ni{t) for the barrier widths b = 4.5, 7.0 nm at T = 10 
K, showing coherent Rabi oscillations as well as subsequent 
scattering-induced relaxation back into the steady state. The 
pump pulse parameters are the same as for the pump-probe 
results (the pump pulse is shown by the dashed line in the 
figures) . 

tering period Tgcatt, and thus charge oscillations cannot 
be observed. Additionally, there are fast oscillations in 
the active region on a time scale of around 300 fs which 
we already addressed in the discussion of the pump-probe 
signals. The oscillations in the injector region can be ex- 
plained by coherences between the injector states where 
E2-E5K5 meV = 850 fs. 

It has already been mentioned that the charge oscil- 
lations are not seen in the electron density evolution 
when only the diagonal elements of the density matrix 
are used in Eq. (PO)) . Thus they constitute a coherent 
effect. This is better exemplified by considering directly 
the subband population evolution caused by the strong 
pump pulse excitation. Wc thus consider the dynamics 
of the WS populations n,;(t) = 2/A^j:n.g(t) under ul- 
trafast optical excitation. It should be stressed here that 
the dynamical calculations are carried out with the full 
equation structure including both the diagonal and non- 
diagonal density-matrix elements. Again, the dynamics 
is considered for the same nonlinear excitation used in 
the pump-probe calculations. The results are shown in 
Fig.im 

During the excitation with the pump pulse, the popu- 
lations of the laser levels undergo Rabi flopping, while the 
other levels follow adiabatically j^^i^^i^°i^^ After the exci- 
tation, the populations return to the steady state due 
to the (incoherent) scattering processes. While for b = 
4.5 nm, the Rabi flopping is restricted mainly to the two 
resonantly excited subbands (l',4), the subband flopping 
for b = 7.0 nm includes both the two laser states (4,1') 
and the injector level (5) since both transitions from the 
lower laser state are energetically very close due to the 
very small tunnel splitting. This weakens the oscillator 
strength of each transition, reducing the effective pulse 
area compared to the full population inversion seen for 
a Itt pump pulse in an ideal two-level system. Still, in 
both cases, the gain transition is saturated, leading to 
absorption at this energy. An oscillating inversion of the 
upper und lower laser states corresponding to the oscilla- 
tions in the pump-probe signals (Fig. or in the spatio- 
temporal electron density evolution on a time scale of a 



few picoseconds (Fig. [TO)) is found neither for b = 4.5 nm 
nor for b = 7.0 nm. As a matter of fact, if the dynamics 
are considered without scattering (not shown), no further 
population dynamics are found after the passage of the 
pulse, while the corresponding pump-probe signals show 
an undamped oscillation in the case of 6 = 4.5 nm. This 
illustrates again that the oscillations in the gain recovery 
arc a coherent effect. 

In summary, it is found that gain oscillations are ob- 
served where the tunnel oscillation period is sufficiently 
less than the scattering period, where the relation 217 > F 
is valid. The tunnel coupling leads to pronounced co- 
herent charge transfer between the injector and the ac- 
tive region as witnessed in the pump-probe signals and 
the spatio-temporally revolved electron density evolu- 
tion. Considerations of quantities which do not directly 
relate to the coherence, such as the dynamics of the WS 
populations, show that the oscillations observed in the 
pump-probe signals constitute a coherent effect which re- 
quires a fully coherent theory to be observed. It should 
be noted here that compared to the current calculations, 
where a stationary localization of charge leads to the drop 
of the current signal, it is here the dynamical transfer of 
charge between different locations in the laser structure 
which constitutes the coherent effect. Thus, the inclu- 
sion of coherence becomes essential in the nonlinear op- 
tical dynamics for the regime where 2Q > F to describe 
the observed coherent effects. Here, the coherence in the 
system drives the gain oscillations, while for 217 < F, the 
scattering destroys the coherent effects. 



IV. CONCLUSION 

We have established a microscopic quantum theory of 
scattering for quantum cascade lasers to treat stationary 
transport as well as nonlinear optical dynamics. A key 
issue is the importance of coherence, i.e. the nondiago- 
nal elements of the density matrix, and the role which it 
plays in different parameter regimes of the tunnel cou- 
pling Q between the two states across the main injection 
barrier and the mean level broadening F of the states. In 
the transport regime, we have compared a rate equation 
(WSH) approach to the current with our microscopic the- 
ory at resonance condition and found that for 2f2 > F, 
the WSH approach is a good approximation. For 217 < F, 
only the full approach shows the physically expected de- 
crease of the current for decreasing tunneling across the 
main injection barrier, whereby the included coherence 
limits the current flowing through the structure via a sta- 
tionary localization of charge in the injection region. On 
the other hand, the consideration of ultrafast nonlinear 
optical pump-probe signals has shown that the inclusion 
of coherence is important to describe the experimentally 
observed coherent effects for 217 > F, where it drives the 
dynamical charge transfer between the injector and the 
active region, resulting in oscillations in the gain recovery. 
For 217 < F, only the typical decrease of the pump-probe 
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signal due to incoherent scattering is observed. It should 
be noted again that, in general, coherences are necessary 
to describe the return to the stationary state after optical 
excitation due to the strong interplay between transport 
and optics in this system. The spatially resolved electron 
density evolution in time shows the density oscillating be- 
tween the injector and the active region for 2^1 > F. The 
consideration of the population dynamics of the levels do 
not show these oscillations, revealing that they are an 
inherent coherent effect resulting from the nondiagonal 
elements of the density matrix. 

Thus, the inclusion of coherence is important in oppo- 
site parameter regimes in stationary transport and non- 
linear optical dynamics: in the former, it becomes impor- 
tant for 217 < r, while in the latter, it allows the observ- 
ability of coherent effects for 2il > T. In this sense, the 
coherence acts in complementary ways: in the transport 
regime, it leads to a limitation of the current which is not 
reproduced by the rate equation approach. In the optical 
regime, the coherence drives the oscillatory modulation 
of the gain recovery which is not destroyed by the scat- 
tering. Thus, it is necessary to consider a fully coherent 



theory in order to describe the combined system of op- 
tics and transport in this nonequilibrium structure in the 
entire parameter range. 
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APPENDIX: SCATTERING EQUATIONS 

1. Electron-LO phonon interaction 

The general form of the scattering equations in Born- 
Markov approximation is given by (for Uq — ^-q) 



ij,k 



(A.l) 



Imn q 



("q + l)'5(e„fc - £Lo) + "q'5(e„fc 



X |(^rtj fnj,k)flm,k+q 
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+ eLo) + nqSie^^j: - e,^^:^- - elo) 



where A™,j(q) = gfjgmn- t^ic calculations, the con- 
siderations are restricted to terms linear in the density 
matrix since in the systems considered here, l/yjl ^ 1 
so that If., rf, rl < 1/ ■ rl- 

The Boltzmann relaxation is determined via 



• (0) 

71 - 

ik 



ik ik 



with the semiclassical Boltzmann in- and out-scattering 
rates given by 



ik 
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2. Ionized doping centers - Impurity scattering 



eLo)l (A-4) by 



The general form of the scattering equations is given 
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4,^ = -f EE Mr^(g)/„,,,5(, 



rnk-\-q-' jn \^ / J nin,k-\-q ^ jk nk-\-q-' 

nm\^/J'in,k ^ nk ink~\-q^ in nm,k-\-q ^ tk nk-\-q' 

I 



(A.5) 



where the average M,^„(g) = (Er,/' ('?)K5m(7'7)>s 
over the statistically distributed doping atoms is intro- 
duced. Since the semiclassical interaction is linear in the 
density matrix, all terms are taken along for this inter- 
action. 

The semiclassical Boltzmann in- and out-scattering 
rates are given by 



pin = ^ y Mll5{e.r - 

ik h ^-^ Ik+q' 



(A.6) 
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